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HIGHLIGHTS 


•  Optimal  control  strategies  for  renewable  hydrogen  production  are  studied. 

•  AID  dynamic  solid  oxide  electrolyser  model  is  coupled  with  a  compressor  model. 

•  The  effects  of  electricity  input  from  intermittent  renewable  energies  are  examined. 

•  Three  different  control  strategies  are  investigated,  with  relative  merits  revealed. 

•  The  energy  consumption  of  the  SOE  system  is  around  35  kWh  per  kg  H2. 
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The  penetration  of  intermittent  renewable  energies  requires  the  development  of  energy  storage  tech¬ 
nologies.  High  temperature  electrolysis  using  solid  oxide  electrolyser  cells  (SOECs)  as  a  potential  energy 
storage  technology,  provides  the  prospect  of  a  cost-effective  and  energy  efficient  route  to  clean  hydrogen 
production.  The  development  of  optimal  control  strategies  when  SOEC  systems  are  coupled  with 
intermittent  renewable  energies  is  discussed.  Hydrogen  production  is  examined  in  relation  to  energy 
consumption.  Control  strategies  considered  include  maximizing  hydrogen  production,  minimizing  SOEC 
energy  consumption  and  minimizing  compressor  energy  consumption.  Optimal  control  trajectories  of 
the  operating  variables  over  a  given  period  of  time  show  feasible  control  for  the  chosen  situations. 
Temperature  control  of  the  SOEC  stack  is  ensured  via  constraints  on  the  overall  temperature  difference 
across  the  cell  and  the  local  temperature  gradient  within  the  SOEC  stack,  to  link  materials  properties 
with  system  performance;  these  constraints  are  successfully  managed.  The  relative  merits  of  the  optimal 
control  strategies  are  analyzed. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Renewable  energy  resources  have  undergone  spectacular 
development  in  recent  years,  to  reduce  our  dependence  on  fossil 
fuels  and  to  mitigate  its  associated  pollutant  gas  emissions  and 
climate  change.  The  existing  worldwide  renewable  power  capacity 
reached  about  1470  GW  in  2012,  up  almost  8.5%  from  2011.  In 
particular,  the  total  global  wind  energy  capacity  and  solar  photo¬ 
voltaic  (PV)  capacity  reached  293  GW  and  100  GW  respectively  in 
2012,  being  three  times  the  installed  wind  capacity  and  ten  times 
the  PV  capacity,  respectively,  that  was  available  just  five  years 
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earlier  [1],  An  electricity  system  based  on  intermittent  renewable 
energy  sources,  such  as  solar,  wind,  hydropower,  geothermal  and 
marine,  gives  rise  to  new  challenges  concerning  the  storage  and 
utilization  of  surplus  energy,  system  operation,  energy  supply 
reliability,  and  future  integration  with  an  electrified  automotive 
sector.  Energy  storage  is  essential  to  help  cope  with  these  chal¬ 
lenges  and  facilitate  a  faster  penetration  of  intermittent  energy 
resources.  Alongside  the  several  proposed  energy  storage  tech¬ 
nologies  such  as  thermal  energy  storage,  battery  and  super¬ 
capacitor  storage,  energy  storage  in  the  form  of  hydrogen  is 
becoming  of  increasing  interest. 

Among  the  many  hydrogen  production  techniques,  water  elec¬ 
trolysis  is  probably  the  cleanest  when  combined  with  a  renewable 
energy  source  to  produce  high  purity  hydrogen.  Although 
commercially  available,  low  temperature  water  electrolysis  based 
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on  alkaline  and  proton  exchange  membrane  (PEM)  electrolysers 
has  not  had  a  significant  penetration,  mainly  due  to  its  high  elec¬ 
tricity  consumption  and  the  challenges  arising  from  cycling  when 
coupled  with  renewables.  High  temperature  (500—1000  °C)  elec¬ 
trolysis  for  hydrogen  production  using  a  solid  oxide  electrolysis  cell 
(SOEC)  appears  to  be  a  promising  technique  because  it  offers 
reduced  electrical  energy  consumption  per  unit  of  hydrogen 
compared  to  low  temperature  electrolysis  due  to  the  favourable 
combination  of  thermodynamics  and  kinetics  at  high  temperatures 
[2,3], 

An  SOEC  can  be  seen  as  the  reverse  of  a  solid  oxide  fuel  cell 
(SOFC),  implying  that  the  development  of  SOEC  technology  can  take 
advantage  of  SOFC  technology  development.  However,  the  under¬ 
lying  electrochemical  mechanisms  in  an  SOEC  and  an  SOFC  are 
likely  to  be  different,  as  is  the  operation  and  application  context.  In 
fact,  experimental  testing  of  SOFC  cells  under  an  SOEC  operating 
mode  has  revealed  degradation  that  is  not  observed  with  the  SOFC 
operating  mode.  [4]  This  has  driven  experimental  research  towards 
understanding  the  microstructure  and  operational  effects,  as  well 
as  developing  new  materials  and  stack  designs  to  deal  with  strong 
reducing/oxidising  environment  at  high  operational  temperature 
[5—13],  Computational  models  have  also  been  developed,  aimed  at 
understanding  SOEC  performance  in  relation  to  operating  condi¬ 
tions  and  material  properties,  and  also  to  assist  the  optimization  of 
SOEC  performance.  These  computational  models  are  mostly  based 
on  physical  and  chemical  descriptions  of  the  SOEC  device  itself. 
These  include:  (1)  the  zero-dimensional  (0D)  steady  state 
isothermal  model  of  a  cell  developed  by  Ni  et  al.  [14]  to  predict  the 
steady-state  electrochemical  behaviour;  (2)  the  one-dimensional 
(ID)  models  developed  by  Grondin  et  al.  [15]  to  predict  the 
steady-state  voltage— current  (V— I)  behaviour  and  the  temperature 
distribution  in  a  cell,  and  Udagawa  et  al.  [16,17]  to  predict  the  dy¬ 
namic  behaviour  of  the  voltage  and  temperature  of  an  SOEC  stack; 
(3)  the  two-dimensional  (2D)  steady-state  model  of  a  single 
repeating  unit  of  a  stack  developed  by  Laurencin  et  al.  [18  to 
predict  V— 1  behaviour  and  distributions  of  temperature,  heat 
fluxes,  local  current  densities  and  gas  concentrations;  and  (4)  the 
three-dimensional  (3D)  computational  fluid  dynamics  (CFD)  model 
of  a  multi-cell  stack  developed  by  Hawkes  et  al.  [19]  to  compute 
spatial-distributed  steady-state  profiles  of  temperature,  electrical 
potential  and  gas  composition. 

The  control  strategy  used  during  operation  is  an  important 
element  which  should  be  taken  into  account  when  designing  and 
defining  the  operational  characteristics  of  the  SOEC  system.  It  is 
even  more  important  when  the  SOEC  system  is  coupled  with 
intermittent  energy  resources  for  hydrogen  production,  due  to  the 
time-varying  nature  of  such  energy  sources.  Although  a  limited 
number  of  system  optimization  studies  have  been  carried  out  to 
investigate  control  strategies  for  SOFC  systems  (which  may  include 
a  fuel  cell  stack,  a  fuel  processing  unit  or  reformer,  air  compressor, 
regulating  pressure  valves,  heat  exchangers  etc.)  [20—23],  to  the 
authors'  knowledge  there  has  not  been  such  a  study  for  SOEC  sys¬ 
tems.  The  only  system-level  modelling  study  for  the  SOEC  system 
was  reported  by  O’Brien  et  al.  [24]  who  presented  a  techno- 
economic  analysis  of  the  SOEC  system  in  the  context  of  its 
coupling  with  nuclear  power  to  explore  the  feasibility  of  different 
system  configurations.  For  such  studies,  lumped-parameter  models 
are  often  used  to  represent  the  behaviour  of  the  devices  involved. 

In  this  paper,  we  present  our  preliminary  results  on  the  iden¬ 
tification  of  optimal  control  strategies  for  the  operation  of  SOEC 
systems  to  offer  efficient  large-scale  operation.  Three  different 
control  objectives  relevant  to  system  operation  are  investigated 
and  discussed.  The  objective  of  this  study  is  to  examine  hydrogen 
production  in  relation  to  energy  consumption,  whilst  maintaining 
control  over  the  SOEC  stack  temperature  distribution,  by  employing 


a  ID  dynamic  SOEC  stack  model  [17]  within  the  overall  system 
model,  to  allow  a  detailed  examination  of  the  performance  of  the 
stack.  The  description  of  the  system  model  and  the  optimal  control 
strategies  are  given  in  the  next  section,  followed  by  results  and 
discussions  in  Section  3.  Conclusions  from  this  study  are  given  in 
Section  4. 

2.  Methodology 

A  complete  SOEC  system  has  to  include  an  energy  source  to 
provide  the  electrical  input,  an  SOEC  stack  as  the  core  part,  and  a 
number  of  auxiliary  devices  (e.g.,  compressor,  reactant  flow  regu¬ 
lator,  heat  exchanger,  steam/hydrogen  separator,  thermal  man¬ 
agement  unit,  power  management  unit).  Precise  control  of  the 
whole  system  is  critical  to  the  viability,  efficiency  and  robustness  of 
the  SOEC  system.  However,  development  of  control  strategies  for 
the  whole  system  is  a  complex  and  computationally  expensive  task, 
as  it  requires  models  for  all  the  devices  involved,  to  explore 
quantitatively  the  interactions  and  trade-offs  between  the  different 
units. 

The  aim  of  this  study  is  to  examine  hydrogen  production  in 
relation  to  energy  consumption.  Two  system  components  domi¬ 
nate  the  energy  consumption  of  the  system,  the  SOEC  stack  and  the 
air  compressor  (when  it  is  used  to  blow  air  for  temperature  control 
[22]).  As  such,  a  simplified  SOEC  system  model  is  explored  here  to 
include  an  SOEC  model  and  a  compressor  model,  as  depicted  in 
Fig.  1.  A  ID  dynamic  SOEC  model  previously  developed  [17]  is 
employed  to  allow  the  detailed  examination  of  the  hydrogen  pro¬ 
duction  associated  with  different  operating  conditions,  and  the 
control  of  the  operation  of  the  SOEC.  Particularly,  temperature 
control  of  the  SOEC  stack  is  only  possible  using  a  dynamic  SOEC 
model  of  at  least  ID,  as  it  allows  the  temperature  gradient  within 
the  cell  to  be  examined,  and  enables  temperature  control  to  avoid 
thermal  excursion  and  cell  failure.  The  electrical  input  is  allowed  to 
vary  over  time  to  represent  the  coupling  of  the  SOEC  system  with 
intermittent  energy  resources. 

Detailed  descriptions  of  the  ID  dynamic  SOEC  model  and  the  air 
compressor  model  are  given  in  this  section,  together  with  expla¬ 
nations  of  the  control  strategies  studied.  The  advantages  and  dis¬ 
advantages  of  each  strategy  are  also  discussed. 

Note  that  in  this  paper,  we  model  a  single  unit  cell  in  the  middle 
of  an  SOEC  stack  to  represent  the  performance  of  the  stack, 
neglecting  end  effects.  Therefore,  air  flow  from  the  compressor  and 
the  electricity  input  to  the  system  are  tailored  for  a  unit  cell,  as  are 
the  hydrogen  production  and  energy  consumption  results  to  be 
reported  later.  This  has  the  benefit  of  having  a  representative  case 
study  which  is  not  limited  by  the  size  of  the  system,  so  that  the 
values  (particularly,  hydrogen  production  and  energy  consump¬ 
tion)  reported  in  this  paper  can  easily  be  scaled  to  systems  of  any 
size,  based  on  the  number  of  the  cells  in  the  SOEC  stack,  assuming 
that  every  cell  in  the  stack  gives  the  same  performance. 


Electricity 


Fig.  1.  Schematic  illustration  of  the  SOEC  system  which  includes  an  SOEC  stack  for 
hydrogen  production  and  an  air  compressor  to  provide  air  for  temperature  control  of 
the  SOEC  stack. 
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2.1.  Description  of  the  ID  dynamic  SOEC  model 

The  ID  dynamic  SOEC  model  developed  is  used  to  simulate  a 
unit  cell  in  the  centre  of  a  sufficiently  large  stack,  so  as  to  represent 
the  performance  of  an  SOEC  stack.  A  schematic  view  of  such  a  unit 
cell  is  shown  in  Fig.  2.  A  planar  SOEC  stack  contains  a  number  of 
repeating  cells  to  achieve  a  sufficient  hydrogen  production  rate. 
Each  individual  cell  consists  of  a  three-layer  solid  structure 
(composed  of  porous  cathode,  electrolyte  and  porous  anode)  and 
an  interconnect  plate.  The  gas  channels  are  embedded  in  the 
interconnect.  Steam  is  introduced  at  the  cathode  side  of  the  solid 
structure  where  it  is  reduced  to  hydrogen,  releasing  oxide  ions  in 
the  process.  The  oxide  ions  then  migrate  through  the  electrolyte  to 
the  anode  where  they  combine  to  form  oxygen  molecules,  releasing 
electrons.  For  modelling  purposes,  the  unit  cell  is  considered  to 
consist  of  four  components,  the  cathode  and  anode  gas  streams,  the 
solid  structure  and  the  interconnect.  All  the  gases  involved  are 
assumed  ideal.  The  cathode  stream  inlet  gas  is  composed  of  10% 
mole  H2  and  90%  mole  H2O.  The  addition  of  hydrogen  in  the 
cathode  stream  mitigates  the  oxidation  of  materials  that  might  be 
induced  by  using  pure  steam.  On  the  anode  side,  air  (modelled  as  a 
mixture  of  21%  mole  02  and  79%  mole  N2)  is  introduced  co- 
currently  with  the  cathode  stream  to  enable  temperature  control 
of  the  stack.  The  stack  is  designed  to  operate  around  1073  K  [3,4]; 
the  temperature  variation  that  occurs  during  stack  operation  is 
controlled  by  introducing  hot  air  (T  =  1073  K)  on  the  anode  side  and 
manipulating  the  air  flow  rate,  as  indicated  in  Section  2.1.4.  Please 
note  that  the  exact  nature  of  the  heat  supply  is  not  considered  in 
this  study;  we  assume  that  waste  heat  from  a  separate  process  is 
used.  The  operating  pressure  is  1.2  x  105  Pa.  The  inlet  flow  rate  of 
the  cathode  stream  is  /c°  =  2  x  10~3  mol  s-1  which  is  equivalent  to 
8920  ml  min-1. 

The  model  is  used  to  simulate  a  cell  with  dimensions  of 
10  cm  x  10  cm  in  the  x  and  y  directions,  with  variations  in  the  y 
direction,  perpendicular  to  the  flow  direction,  being  ignored.  The  z 
direction  gives  the  heights  of  the  cell  components  and  also  the 
height  of  the  whole  stack.  A  detailed  description  of  the  cell  ge¬ 
ometry  and  material  properties  can  be  found  in  Appendix  [25],  The 
equations  are  developed  and  solved  along  the  x  direction.  The  re¬ 
sults  to  be  presented  are  normalized  along  the  cell  length,  with 
x  =  0  indicating  the  inlet  and  x  =  1  indicating  the  outlet. 

The  model  developed  consists  of  an  electrochemical  model, 
mass  balances  for  the  gas  streams,  and  energy  balances  for  the 
cathode  and  anode  streams,  solid  structure  and  interconnect.  The 
equations  are  given  in  Tables  1—3,  and  the  nomenclature  can  be 
found  at  the  end  of  this  paper. 
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Fig.  2.  Schematic  view  of  a  unit  cell  of  a  planar  SOEC  stack. 


X=0 


o2/n2  — -  o2/n2 


Interconnect 

H,  /  H2C)  Cathode  str  eam  -*■  ft,  /  H2Q 


Solid  structure 


G;  /  N2  — Anode  stream  — 02  /  N2 

Interconnect 


2.1.1.  Electrochemical  model 

An  electrochemical  model  relates  variables  such  as  gas  species 
concentrations,  cell  component  temperatures  and  average  current 
density  to  the  electrical  potential  of  the  cell,  which  can  then  be  used 
to  calculate  the  electrical  energy  consumption.  The  cell  potential 
corresponds  to  the  sum  of  the  reversible  potential  and  all  the 
irreversible  losses  that  occur  as  the  electrical  current  passes 
through  the  cell.  Irreversible  losses  consist  of  cathode  and  anode 
activation  overpotentials,  cathode  and  anode  concentration  over¬ 
potentials  and  ohmic  losses.  These  irreversible  losses  are  partly 
responsible  for  the  heat  produced  within  the  cell.  Mathematical 
descriptions  of  the  electrochemical  model  are  given  in  Table  1. 


2.1.2.  Mass  balances 

The  mass  balances  represent  the  gas  flow  behaviour  inside  the 
cathode  and  anode.  The  mass  balance  equations  are  based  on  the 
mass  conservation  principle,  and  the  thermodynamic  properties  of 
gas  mixtures.  The  composition  of  the  cathode  stream  changes  along 
the  cell,  and  is  taken  into  account  in  the  model  using  the  cathode 
stream  mass  balances  and  boundary  conditions  at  the  cell  inlet, 
which  allows  the  prediction  of  hydrogen  and  water  concentrations 
at  each  position  along  the  channel.  For  the  anode  stream,  mass 
balances  and  boundary  conditions  at  the  inlet  predict  the  O2  and  N2 
concentrations  at  each  position  along  the  cell.  Mathematical  de¬ 
scriptions  of  the  mass  balances  are  given  in  Table  2. 


2.1.3.  Energy  balances 

Energy  balances  are  a  key  driver  to  predict  the  temperature 
distribution  along  the  cathode  and  anode  streams,  solid  structure 
and  interconnect.  Heat  exchange  between  the  gas  streams  and 
walls  (solid  part  of  the  cells)  is  fully  described  by  convection.  The 
thermal  fluxes  along  the  solid  parts  of  the  cell  can  be  modelled 
using  Fourier's  law  of  heat  conduction  while  taking  into  account  the 
radiative  heat  exchange  between  walls  and  interconnect.  Mathe¬ 
matical  descriptions  of  the  energy  balances  are  given  in  Table  3. 


2.1.4.  Temperature  control 

Temperature  control  is  provided  for  the  SOEC  stack  through 
manipulation  of  the  air  ratio  on  the  anode  side.  The  air  ratio  pro¬ 
vides  a  measure  of  the  inlet  flow  rate  of  air  in  relation  to  the  rate  of 
reaction.  It  is  defined  as  the  ratio  between  the  molar  flow  rate  of  O2 
in  the  inlet  air  flow  to  the  number  of  moles  of  O2  produced  in  the 
unit  cell,  per  unit  time.  The  mathematical  expression  for  the  air 
ratio,  is 


2F/AV)yO2(0.^ 

j(t)Lwr02 


(1) 


As  indicated  by  Eq.  (1 ),  the  inlet  flow  rate  of  the  anode  stream 
is  co-determined  by  the  air  ratio  and  operating  current  density,  and 
varies  with  operating  conditions. 


2.2.  The  compressor  model 

The  compressor  model  in  this  study  is  a  simple  steady-state 
model  which  links  the  air  compressor  power  to  the  air  stream 
flow  rate,  ma ir,  and  temperature,  T,  the  caloric  properties  of  air  (i.e. 
cp  and  y),  and  the  inlet  and  outlet  of  the  pressure  through  the 
compressor. 

The  parasitic  power  required  by  the  air  compressor,  which 
supplies  the  total  amount  of  air  needed  to  meet  the  stack  cooling 
requirement,  is  modelled  as  follows  [22]: 


Table  1 

The  dynamic  SOEC  model:  electrochemical  model. 
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Cell  working  potential:  U(t)  =  Urev(X,  t )  +  770hm(X,  0  +  *?conc, cathode  (*,  t)  +  *7conc, anode  (*,  0  +  ^act, cathode (*>  0  +  *7act, anode (*,  t),  te  [0,  t/] 


Reversible  potential :  Urev(x,  t)  =  U°(x ,  t) 


XTs(x,t)  ln  f  (x,t)(p0?  (t)xlQ— 


Ch2o(x,0 


,  te  [0,  t/] 


Concentration 


^cathode  electrolyte  (*if)  ffanode 

overpotential:  T7conC|Cathode (*> 0  =  In  fe ((* :  ^conc, anode  (*>  0  = 


C™(x,r)Ts(x,t) 

,t)TA(x,t) 


C™(x,  t)  =  Ch2  (x,  t)  +a£^Lrj(x,  t):  C™  (x,  t)  =  CH2o(x,  t)  -  t) 

C™(x.  t)  =  C0i  (*,  t)  +  CN2  (x,  t)  -  CN2  (x,  t)exp (x.q])  ■  te I0'  V 


Activation  overpotential: 


j(*>  0  —  Jo, cathode  f) 


exp f 2(/.~tt)F 77  (x  5B°!x,t! exp f  ~2“F,T7  .  .  (x  t) 

Ch2  (x,t)  ex”  l  5RTs(x,t)  ^act, cathode  1 1  J  CH*o(x,0  ex”  l  MTS (x,t)  4 act, cathode ^  L) 


,  te[0,t/] 


j(^>  Jo, anode  (*^5  0 


eXP  (#rs(x,0  ^act, anode  (x>  0^  ~  C?2  (x,t)  eXP  ^act, anode  (x> 


Jo, electrode  (x,t)  =  kdedrode  exP  [n^jf]  ’  electrodes  {cathode,  anode},  tefO,^] 


,  telO.tf] 


Table  2 

The  dynamic  SOEC  model:  mass  balances. 


factor  of  air  pa jr  and  the  molar  volume  ( RT/P )  for  the  purpose  of  unit 
conversion. 


Cathode  stream  mass  balance 

Jj[Q(x,t)]  =  -uCi|[Ci(x,t)]+^^(x,t),  ie{H2,H20},  ts[0 ,tf] 
Cj(0,0)  =  C°,  ie{H2,H20} 

R(x,t)  =  $$,  ts  [0,  tf] 

Anode  stream  mass  balance 

|[Q(x,t)]  =  -uAi[Ci(x,t)]  +  ^./1R(x,t),  ie{02,N2},  te[0,tjr] 
Q(0,0)  =  C,P,  is{02,N2},  te[0,tf] 


mair(f)  =  Pair(I',P)j?(t)T- 


2.3.  Optimal  control  strategies 
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Pcp(t) 


mair(t) 


CP,3lJ(t) 

^EM^lCP 
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The  inlet  air  pressure  p jn  is  equal  to  the  ambient  pressure  of 
1.0  x  105  Pa.  The  outlet  air  pressure  pout  is  equal  to  the  operating 
pressure  of  the  SOEC  stack,  which  in  this  case  is  1.2  x  105  Pa.  This 
gives  a  compression  ratio  p0ut/Pm  =  1.2. 

The  air  mass  flow  rate  mair  in  Eq.  (2)  is  linked  to  the  inlet  flow 
rate  of  the  anode  stream  jj(t)  in  the  SOEC  by  adding  the  density 


The  development  of  suitable  optimal  control  strategies  for  the 
SOEC  system  relies  on  an  understanding  of  the  characteristics  of 
both  the  SOEC  system  and  the  coupled  electricity  supply.  Two 
scenarios  are  considered  in  this  study. 

A.  The  heating  and  cooling  of  an  SOEC  stack,  during  which  the 
current  density  is  a  manipulated  variable. 

B.  The  SOEC  system  is  integrated  with  an  intermittent  source  of 
energy,  during  which  the  electrical  energy  input  is  treated  as  a 
disturbance. 

In  Scenario  A,  the  optimal  control  unit  is  used  to  set,  within  a 
given  range  of  current  densities,  the  current  density  at  which  the 


Table  3 

The  dynamic  SOEC  model:  energy  balances. 


Cathode  stream 

»  [Tc(x,t)l  =  -uCi^  [Tc(x,  t)]  +  [TS(X,  t)  -  Tc(x,t)]  +-Jkjj_.irI(x,t)  -  Tc(x,t)],  t<=[0,  tf] 

Tc( 0, 0 )  =  T°,kc  =  Nuc£-c,  dh, c  =  t<=  [0,  tf] 

Anode  stream 

„  JL  \Ta  (X  til  4  ,  ...  .  , 

PACp,A"A  1  0  '  n '  /J  PACp.Atl, 

Ta(0, 0)  =  Ti(Q’0)+T'(0’0).  kA  =  Nua^,  dhA  =  ^,  te [0,  tf] 

Solid  structure 

m  Ps(*.  01  =  sti  w  Ps(*>  f)I — 81 


m  Wa(x,  t)]  =  -uAJj  [Ta(x,  t)]  +  — ^  [Ts(x,  t)  -  rA(x,  t)]  +  — [T|(x,  t)  -  TA(x,  t)],  te  [0,  tf] 


*cP,s*s  Ps(x.  0  -  Tc(x,  t)]  - [Ts(x,  t)  -  TA(x,  t)]  -  ^ 
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-AH(x,  t)R(x,  t )  +  j(x,  t)U(t )],  te  [0,  tf\ 


fx[Ts( 0,0)]  =  0,  ^[Ts(l,0)]  =  0,  te[0,t/] 

Interconnect 
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Table  4 

Mathematical  descriptions  of  optimal  control  strategies. 


Representative 

scenarios 

Optimization  objective 

Constraints 

Control  variables 

A.  Heating  and 

(a)  Maximizing  hydrogen  production 

Path  constraints 

•  Time  horizon  tf 

cooling 

1  Acceptable  temperature  gradients 

•  Air  ratio  0.4  <  \{/  <  14 

Th2, total  -  /  YH!(t)dt  Yjt(t)  =/c(t)CH2(l,t)  -/c°(t)CH2(0,t) 

across  entire  SOEC 

•  Current  density 

(b)  Minimizing  SOEC  energy  consumption 

ft=tf 

Esoec  =  /  PsoEc(0dt  Psoec(0  =  U(tjj(t) 

Jt= 0 

|ATc(t)|  =  |Tc(l,t)-rc(0,t)|  <  100/C.  te[0,t/] 

2  Acceptable  local  temperature  gradients 
|3Ts(x,t)/0x|  <  10/C,  te[0,t/] 

3  Limiting  power  to  compressor 
o  <  PcM  <  PcPmax.  fe[0,t/] 

End-point  constraints 

4  Returning  to  steady-state 

9ATc{t)/at  =  0  at  t  =  tf 

1000  <]<  7000  Am  2 

B.  Electricity 

(c)  Minimizing  compressor  energy  consumption 

Path  constraints 

•  Time  horizon  tf 

disturbance 

1  Acceptable  temperature  gradients 

•  Air  ratio  0.4  <  \f/  <  14 

Ecp  =  /  Fcp(t)dt 

|ATc(t)|  =  |Tc(l,t)-rc(0,t)|  <  100IC,  te[0,tf] 

•  Total  input  power 

2  Acceptable  local  temperature  gradients 
|0Ts(x,t)/0x|  <  10/C,  te[ 0,  tf] 

End-point  constraints 

3  Returning  to  steady-state 
dATc{t)/dt  =  0  at  t  = 

20  <  Ptotal  <  100  W 

SOEC  stack  should  be  operating  in  order  to  achieve  the  control 
objectives.  The  electrical  energy  demand  in  this  scenario  is  thus 
determined  by  how  the  SOEC  stack  is  operated,  in  Scenario  B, 
however,  the  electrical  energy  availability  is  pre-determined  as  an 
input  to  the  SOEC  stack;  the  optimal  control  unit  reacts  accordingly 
under  this  electricity-input  disturbance  to  achieve  the  control 
objectives. 

2.3.1.  Control  objectives 

Control  objectives,  which  reflect  the  purposes  and  functions  of  a 
system,  can  be  chosen  depending  on  the  specific  operating  context. 
Three  optimal  control  situations  are  considered  with  different 
control  objectives:  (a)  maximizing  hydrogen  production,  (b)  mini¬ 
mizing  the  electrical  energy  consumption  from  the  SOEC  stack,  and 
(c)  minimizing  the  electrical  energy  consumption  from  the  air 
compressor.  The  first  two  optimal  control  situations  are  studied 
under  Scenario  A,  whilst  the  third  optimal  control  situation  is 
studied  under  Scenario  B.  For  all  the  control  situations  considered, 
certain  constraints  have  to  be  satisfied  in  order  to  keep  the  system 
working.  This  is  further  discussed  in  the  next  sections. 

2.3. 7.1.  Scenario  A:  heating  and  cooling 

23.1.1.1.  Maximizing  hydrogen  production.  This  control  strategy 
may  arise  when  there  is  abundant  electricity  available  from  inter¬ 
mittent  energy  sources,  for  instance  in  the  case  when  the  electricity 
generated  exceeds  the  demand.  In  this  situation,  maximizing 
hydrogen  production  is  desirable  so  that  the  surplus  electricity  can 
be  stored  in  the  form  of  hydrogen.  This  control  strategy  may  also  be 
suited  when  there  is  a  high  hydrogen  demand  in  the  supply  chain. 
The  amount  of  hydrogen  produced  is  directly  dependent  on  the 
current  density  at  which  the  SOEC  is  operating.  As  an  example,  here 
we  set  lower  and  upper  bounds  on  the  operational  current  density 
range  (see  Table  4),  and  let  the  optimal  controller  determine  how 
the  SOEC  should  be  operating,  in  accordance  with  Scenario  A.  The 
objective  function  for  maximizing  the  hydrogen  production  of  the 
process  over  a  given  time  horizon  is  given  in  Table  4.  The  hydrogen 
production  rate,  YH2,  is  calculated  as  the  molar  flow  rate  of 
hydrogen  at  the  outlet  of  the  cathode  stream  minus  the  molar  flow 
rate  of  hydrogen  at  the  inlet  of  the  cathode  stream. 

23.1.1.2.  Minimizing  the  energy  consumption  of  the  SOEC  stack. 
When  the  available  electricity  for  hydrogen  production  is  limited  or 
the  demand  on  hydrogen  production  is  low,  it  is  desirable  to 
operate  the  SOEC  system  such  that  the  electrical  energy 


consumption  is  minimized  to  reduce  the  operating  cost.  As  the 
SOEC  stack  is  one  of  the  main  electrical  energy  consumers  in  the 
system,  it  is  therefore  necessary  to  minimize  the  electrical  energy 
consumption  of  the  SOEC  stack.  The  objective  function  for  mini¬ 
mizing  the  electrical  energy  consumption  of  the  SOEC  stack  within 
a  given  time  horizon  is  given  in  Table  4.  The  power  consumed  by 
the  SOEC  per  unit  time  Psoec  is  calculated  as  the  product  of  the  cell 
voltage  and  the  average  current  density  of  the  SOEC  stack.  To 
compare  with  case  (a)  of  maximizing  hydrogen  production,  this 
optimal  control  situation  is  studied  with  the  same  control  variables 
as  (a),  where  the  optimal  controller  determines  how  the  SOEC 
should  be  operating  given  an  operational  current  density  range 
with  fixed  lower  and  upper  bounds. 

23.1.2.  Scenario  B:  electricity  disturbance 

23.1.2.1.  Minimizing  electrical  energy  consumption  from  air 
compressor.  Here  we  consider  control  under  Scenario  B  where  the 
electrical  energy  is  a  disturbance  to  the  system.  In  this  scenario,  the 
possibility  of  both  maximizing  hydrogen  production  and  mini¬ 
mizing  electrical  energy  consumption  of  the  SOEC  system  at  the 
same  time  is  explored.  To  achieve  both  objectives  simultaneously, 
the  control  objective  is  set  to  minimize  the  electrical  energy  con¬ 
sumption  from  the  air  compressor.  For  a  given  electricity  input,  by 
minimizing  the  electrical  energy  consumption  of  the  air 
compressor,  more  electrical  energy  is  then  available  for  the  SOEC 
stack  which  will,  in  turn,  maximize  the  amount  of  hydrogen  pro¬ 
duced.  The  objective  function  for  minimizing  energy  consumption 
from  the  air  compressor  is  given  in  Table  4. 

23.2.  Constraints 

In  optimal  control  problems  it  is  vital  to  identify  the  constraints 
for  a  given  process,  as  constraints  determine  the  limit  of  the 
feasibility  and  the  flexibility  of  the  process.  In  the  optimal  control 
studies  presented  in  this  paper,  two  kinds  of  constraints  are 
imposed,  namely  path  constraints,  and  end-point  constraints.  Path 
constraints  have  to  be  satisfied  at  all  times  during  operation,  whilst 
end-point  constraints  have  to  be  satisfied  only  at  the  end  of  the 
time  horizon.  The  mathematical  forms  of  the  constraints  are  given 
in  Table  4. 

For  all  the  control  situations,  it  is  important  to  ensure  that  the 
temperature  profile  within  the  SOEC  stack  satisfies  the  constraints 
on  the  temperature  gradients  at  all  times  to  avoid  excessive  ther¬ 
mal  excursion  of  the  ceramic  materials.  In  previous  studies,  a 
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maximum  allowable  local  thermal  gradient  along  the  cell  (partic¬ 
ularly  in  the  solid  structure  which  consists  of  cathode,  electrolyte 
and  anode)  has  been  estimated  to  be  10  K  citT1,  and  a  maximum 
allowable  overall  temperature  difference  along  the  cathode  stream 
for  a  cell  of  10  cm  to  be  100  K  [26].  Two  temperature  constraints  are 
thus  implemented:  the  absolute  temperature  difference  | ATC| 
across  the  10  cm  cell  must  be  less  than  or  equal  to  100  K  (constraint 
1  in  Table  4);  and  the  absolute  local  temperature  gradient  |3Ts/3x| 
must  be  less  than  or  equal  to  10  K  cm-1  (constraint  2  in  Table  4).  In 
addition,  an  end-point  constraint  is  imposed  on  the  derivative  of 
the  overall  temperature  difference  to  ensure  that  the  system  rea¬ 
ches  equilibrium  (i.e.  no  fluctuation)  at  the  end  of  the  operation 
(constraint  4  in  Scenario  A  and  3  in  Scenario  B,  Table  4). 

Furthermore,  in  Scenario  A,  considering  that  the  air  compressor 
is  one  of  the  main  electrical  energy  consumers,  it  is  also  desirable  to 
limit  its  electrical  energy  consumption  whilst  fulfilling  the  tem¬ 
perature  control  requirements  of  the  SOEC  stack.  A  path  constraint 
is  imposed  to  provide  an  upper  bound  over  the  electrical  energy 
consumption  by  the  air  compressor.  Its  formulation  is  shown  as 
constraint  3  in  Scenario  A.  Its  effects  will  be  shown  and  discussed  in 
the  results  section.  Pmax  is  adjusted  according  to  the  system  oper¬ 
ating  conditions,  in  relation  to  current  density  and  air  flow  rate,  as 
indicated  by  Eqs.  (1)  and  (2).  This  helps  us  explore  trade-offs  with 
the  air  compressor  power  control.  If  any  of  the  constraints  is  not 
satisfied  during  operation,  the  process  is  infeasible. 


Fig.  3.  A  schematic  of  the  control  variable  profiles  that  can  be  obtained  with  control 
vector  parameterization,  for  piecewise  constant  and  piecewise  linear 
parameterizations. 

Both  the  durations  of  the  time  intervals  and  the  values  of  the 
control  parameters  are  treated  as  variables. 

The  optimal  control  profiles  reported  in  this  paper  are  obtained 
with  piecewise  constant  controller;  a  piecewise  linear  controller 
was  also  investigated,  but  it  does  not  bring  benefits  in  terms  of 
giving  a  better  solution. 


2.3.3.  Control  variables 

In  all  the  optimal  control  situations  considered,  the  control  time 
horizon,  t/,  is  set  to  be  3000  s.  The  air  ratio,  defined  as  the  ratio  of 
the  inlet  flow  rate  of  air  on  the  anode  side  to  the  gas  inlet  flow  rate 
on  the  cathode  side,  is  used  as  a  control  variable  to  limit  temper¬ 
ature  variation.  The  lower  and  upper  bound  for  the  air  ratio  are 
selected  to  be  0.4  and  14  respectively,  considering  system  con¬ 
straints  [4],  For  control  situations  (a)  and  (b)  where  we  have  control 
on  the  operational  current  density  of  the  SOEC  stack,  the  operating 
regime  for  current  density  is  between  1000  and  7000  A  nrT2.  For 
control  situation  (c),  where  the  electrical  energy  is  an  input 
disturbance,  various  configurations  of  the  electrical  power  input 
are  selected,  allowing  the  electrical  power  to  vary  between  20  and 
100  W  (for  input  to  a  cell)  over  a  time  horizon  of  3000  s.  The 
mathematical  descriptions  of  the  control  variables  are  given  in 
Table  4. 


3.  Results  and  discussions 

The  optimal  control  problems,  in  which  the  system  of  partial 
differential  and  algebraic  equations  for  the  ID  dynamic  SOEC  stack 
model  the  air  compressor  model  are  embedded,  are  solved  under 
all  scenarios  via  gPROMS  Model  Builder  3.4. 

Before  performing  optimal  control  calculations,  the  steady  state 
performance  of  the  SOEC  is  investigated  to  gain  a  better  under¬ 
standing  of  the  SOEC  performance.  Fig.  4  shows  the  operational  cell 
voltage  and  the  open  circuit  voltage  (OCV)  over  a  wide  range  of 
operating  current  densities  at  1073  K.  When  current  density  in¬ 
creases,  the  cell  voltage  also  increases,  but  remains  below  the 
thermal-neutral  voltage  of  1.30  V  at  1073  K  for  the  entire  range  of 
current  densities  investigated,  featuring  endothermic  operation. 
The  OCV  remains  the  same  over  the  entire  range  of  operating 
current  densities.  When  operating  at  1000,  7000  and  10,000  A  m  2, 
a  voltage  difference  (between  the  operational  voltage  and  the  OCV) 


2.3.4.  Control  vector  parameterization 

Control  vector  parameterization  [27  is  used  to  formulate  the 
optimal  control  problem.  In  this  approach,  control  variables  are 
allowed  to  achieve  solutions  of  a  pre-specified  functional  form,  e.g. 
piecewise  constant  and  piecewise  linear  controls  with  continuity  at 
the  interval  boundaries.  The  time  horizon  is  discretized  into  non- 
uniform  time  intervals.  Two  types  of  parameterization  are  consid¬ 
ered:  piecewise-constant  controls,  for  which  constant  control 
parameter  values  are  assumed  over  each  interval,  and  piecewise- 
linear  controls,  for  which  control  parameter  values  with  constant 
positive  or  negative  gradients  are  assumed  over  each  interval,  as 
illustrated  in  Fig.  3.  More  complex  forms  can  also  be  considered. 

The  benefit  of  this  approach  is  that  the  infinite-dimensional 
optimal  control  problem  is  reduced  to  an  optimization  problem 
in  a  few  variables,  e.g.  in  the  case  of  piecewise  constant  control,  a 
control  variable  u(t)  is  described  by: 


(  uu  0  <  t<  tj 
I  u2,  0  <t<t2 

{  uf,  0  <t<tf 


(4) 


Fig.  4.  Operational  cell  voltage  and  open  circuit  voltage  (OCV)  at  a  range  of  current 
densities. 
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of 41, 255,  and  340  mV  is  found  respectively,  indicating  the  increase 
of  irreversible  losses  within  the  cell  as  current  density  is  increased. 
The  area  specific  resistance  (ASR)  calculated  is  0.44  Q  cm2  when 
operating  at  1000  A  itT2,  and  0.36  0  cm2  when  operating  at 
10,000  A  nr2,  showing  a  gradual  decrease  of  ASR  with  increasing 
current  density.  For  Scenario  A,  where  we  have  control  over  the 
average  current  density  j  between  1000  and  7000  A  uT2,  opera¬ 
tions  at  1000  and  7000  A  nT2  are  examined  to  understand  the 
bounds  on  the  air  ratio  that  can  be  used  to  control  the  system  so 
that  the  temperature  constraints  are  met.  Fig.  5  shows  the  cathode 
stream  temperature  distribution  along  the  cell.  As  can  be  seen,  the 
cathode  stream  temperature  decreases  along  the  cell  from  the  inlet 
(normalized  position  x  =  0)  to  the  outlet  (x  =  1,  corresponding  to  a 
length  of  10  cm)  for  all  the  operating  conditions  investigated, 
featuring  endothermic  operation.  When  operating  at  1000  A  m~2 
with  a  minimum  air  ratio  of  0.4,  the  temperature  gradient  between 
the  inlet  and  the  outlet  is  164  K,  which  exceeds  the  maximum 
temperature  gradient  of  100  K  allowed  for  a  10  cm  cell,  and  thus 
puts  the  cell  in  danger  of  breakdown.  A  minimum  air  ratio  of  2.71  is 
required  to  reduce  the  temperature  gradient,  in  this  case,  to  99.3  K. 
The  operation  with  a  maximum  air  ratio  of  14  is  also  examined, 
showing  that  the  temperature  gradient  is  maintained  at  48  K. 
When  operating  at  a  current  density  of  7000  A  itT2,  the  cell  fea¬ 
tures  a  temperature  gradient  of  48  K  along  the  cell  with  a  minimum 
air  ratio  of  0.4,  and  a  temperature  gradient  of  17  K  with  a  maximum 
air  ratio  of  14.  Based  on  the  monotonic  behaviour  of  the  tempera¬ 
ture  gradient  with  the  air  ratio,  this  means  that  any  air  ratio  value  in 
the  control  range  of  0.4—14  is  sufficient  for  the  temperature  control 
to  meet  the  temperature  constraint  at  higher  current  densities. 
From  Fig.  4,  the  initial  guess  for  the  optimization  studies  are  set  to 
be:  air  ratio  of  2.71  when  operating  at  current  density  of 
1000  A  m~2,  air  ratio  of  0.4  when  operating  at  current  density  of 
7000  A  nT2. 

3.3.  Maximizing  hydrogen  production 

Fig.  6  shows  the  optimal  control  profiles  in  the  case  of  maxi¬ 
mizing  hydrogen  production.  To  achieve  maximum  hydrogen  pro¬ 
duction  within  the  control  period  of  3000  s,  current  density 
increases  from  1000  to  7000  A  nT2  within  5  s  at  the  start  of  the 
control  period  and  is  maintained  at  7000  A  nT2  until  the  end  of  the 


Fig.  5.  Cathode  stream  temperature  distribution  along  the  cell.  The  dash-dotted  line 
indicates  operation  at  1000  A  in  2  with  an  air  ratio  of  0.4.  The  dotted  line  indicates 
operation  at  1000  A  m  2  with  an  air  ratio  of  2.71.  The  dashed  line  indicates  operating 
at  1000  A  m  2  with  an  air  ratio  of  14.  The  thick  solid  line  indicates  operating  at 
7000  A  m  2  with  an  air  ratio  of  0.4.  The  thin  solid  line  indicates  operating  at 
7000  A  m  2  with  an  air  ratio  of  14. 


control  period,  as  shown  in  Fig.  6(a).  In  response  to  the  change  in 
the  operating  current  density,  the  cumulative  amount  of  hydrogen 
produced,  Yh2, total >  increases  linearly  with  time  and  reaches  a 
maximum  of  1.09  mol. 

For  operations  at  current  density  of  7000  A  n~r2,  the  imple¬ 
mented  upper  limit  of  the  compressor  power,  Pcpmax.  is  set  to 
different  values  in  the  range  of  3.4—120  W  to  examine  the  effect  of 
temperature  control.  These  values  are  calculated  based  on  the 
operating  current  density)  and  air  ratio  i j/  using  Eqs.  (1)  and  (2); 
Pc Pmax  =  3.4  W  corresponds  to  j  =  7000  A  m~2  and  ^  =  0.4,  while 
Pcpmax  =  120  W  corresponds  to  j  =  7000  A  m~2  and  xf  =  14.  Fig.  6(b) 
shows  how  the  air  ratio  responds  to  the  compressor  power  upper 
bound:  it  rises  or  falls  at  t  =  1  s  from  the  initial  value  of  2.71  to  a 
constant  value  corresponding  to  the  maximum  compressor  power. 
When  Pcpmax  <3.4  W,  it  is  not  feasible  to  control  the  SOEC  as  it 
violates  the  lower  bound  of  the  air  ratio  value  of  0.4;  when  Pcpmax 
>120  W,  the  optimal  control  problem  is  always  feasible  and  the  air 
ratio  always  stays  at  the  upper  bound  of  the  air  ratio  of  14.  Fig.  6(c) 
shows  that  the  temperature  difference  between  the  outlet  and  inlet 
of  the  cell  always  satisfies  the  constraint  |ATC|  <  100  K  and  reaches 
steady  state  towards  the  end  of  the  control  period  (t  =  3000  s).  The 
absolute  value  of  the  maximum  local  temperature  gradient  in  the 
solid  structure,  |9Tj(x,  t)/9x|  is  found  to  be  under  the  upper  limit  of 
10  K  for  all  the  Pcpmax  values  studied,  as  shown  in  Fig.  6(d).  It  is  also 
obvious  that  with  a  small  compressor  power  limit,  the  cell  operates 
with  a  larger  temperature  difference  and  the  temperature  takes 
longer  to  stabilize,  and  vice  versa.  For  example,  with 
Pcpmax  =  120  W,  a  high  air  ratio  of  14  is  used  so  that  a  comparatively 
large  air  flow  rate  is  introduced  to  the  anode  stream  of  the  SOEC 
stack,  bringing  the  initial  temperature  difference  of  ATc  =  -99.3  I< 
to  a  final  temperature  difference  of  A Tc  =  —16.8  K  within  250  s; 
with  Pcpmax  =  4  W,  a  low  air  ratio  of  0.4  is  introduced,  bringing  the 
initial  temperature  difference  of  ATc  =  -99.3  K  to  a  final  temper¬ 
ature  difference  of  A  Tc  =  -46  K  within  1240  s. 

3.2.  Minimizing  the  SOEC  energy  consumption 

Fig.  7  shows  the  optimal  control  profiles  in  the  case  of  mini¬ 
mizing  the  SOEC  energy  consumption.  As  shown  in  Fig.  7(a),  cur¬ 
rent  density  decreases  from  7000  to  1000  A  m~2  at  the  start  of  the 
time  horizon  and  remains  at  1000  A  m  2  till  the  end  of  the  control 
period  (t  =  3000  s).  In  response  to  this,  the  cumulative  energy 
consumption  of  the  SOEC  stack  increases  linearly  with  time  and 
reaches  7.6  Wh  at  t  =  3000  s. 

The  upper  limit  of  the  compressor  power,  Pcpmax.  is  set  to 
different  values  in  the  range  of  3—16  W,  corresponding  to  the  air 
ratio  range  of  0.4—14,  as  calculated  from  Eqs.  (1 )  and  (2).  In  Fig.  7(b), 
one  can  see  how  the  air  ratio  responds  to  the  compressor  power 
bounds:  it  rises  at  t  =  1  s  from  the  initial  value  of  0.4  to  the 
controlled  values  based  on  the  compressor  power  bound,  as  current 
density  changes  from  7000  to  1000  A  nrr2.  For  the  range  of  air  ratio 
values  considered,  the  temperature  difference  between  the  outlet 
and  inlet  of  a  cell  is  controlled  within  the  constraint  |ATC|  <  100  I< 
and  reaches  steady  state  towards  the  end  of  the  control  period 
(t  =  3000  s),  as  shown  in  Fig.  7(c).  The  absolute  value  of  the 
maximum  local  temperature  gradient  in  the  solid  structure, 
|9Ts/9x|,  is  observed  to  be  under  the  10  1<  bound  for  all  the  Pcpmax 
control  cases  as  shown  in  Fig.  7(d).  It  is  also  obvious  that  with  a 
small  compressor  power  limit,  a  larger  temperature  difference  is 
observed,  and  vice  versa.  For  example,  with  Pcpmax  =  16  W,  a  high 
air  ratio  of  14  is  introduced  to  the  cathode  stream  of  the  SOEC  stack, 
which  enables  only  a  slight  change  in  temperature  and  keeps 
A  Tc  — 46  K  through  the  entire  control  period;  with  Pcpmax  =  3  W, 
an  air  ratio  of  2.71  is  introduced,  bringing  the  temperature  differ¬ 
ence  of  ATc  =  -46  K  to  a  constant  value  of  -99.5  K  within  2500  s. 
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Fig.  6.  Optimal  profiles  of  (a)  current  density  and  hydrogen  production,  (b)  air  ratio,  (c)  cathode  stream  temperature  difference  and  (d)  maximum  local  temperature  gradient  in  the 
solid  structure,  for  the  case  of  maximizing  hydrogen  production. 


3.3.  The  effect  of  compressor  power  on  hydrogen  production  and 
energy  consumption 

The  results  in  the  above  two  control  cases  (i.e.  maximizing 
hydrogen  production  and  minimizing  the  SOEC  energy  consump¬ 
tion)  show  that  the  control  of  the  SOEC  system  illustrated  in  Fig.  1  is 
feasible  as  long  as  the  compressor  power  limit  is  maintained  within 
a  reasonable  range,  under  which  the  constraints  become  active;  the 
impact  of  imposing  a  compressor  power  limit  on  the  system  affects 
the  air  flow  rate,  and  consequently  the  temperature  gradient  A Tc 
and  the  rate  of  the  temperature  gradient  change  with  time  dATc/Qt. 
When  the  power  limit  for  the  compressor  is  sufficiently  large,  the 
constraints  become  inactive.  If  a  sufficiently  large  power  limit  is  set 
on  the  compressor,  the  system  will  automatically  choose  to  work 
under  the  biggest  allowable  air  flow  rate,  bringing  ATc  and  c)ATc/3t 
down  to  very  low  values  which  benefit  the  materials  of  the  SOEC 
stack. 

However,  the  power  limit  imposed  on  the  compressor  needs  to 
be  set  considering  other  factors  such  as  the  system  energy  con¬ 
sumption.  The  effects  of  the  compressor  power  limit  on  both  en¬ 
ergy  consumption  and  hydrogen  production  of  the  SOEC  system  are 
further  explored.  Fig.  8(a)  and  (b)  shows  the  total  hydrogen  pro¬ 
duction  and  energy  consumption  over  the  entire  time  horizon  of 
3000  s,  as  a  function  of  the  maximum  compressor  power,  Pcpmax. 
for  both  maximizing  hydrogen  production  and  minimizing  the 
SOEC  energy  consumption.  As  can  be  seen  from  Fig.  8(a)  and  (b),  the 
total  amount  of  hydrogen  produced  does  not  change  with  the 
variation  in  compressor  power,  as  expected.  According  to  Faraday's 


law  of  electrolysis,  the  amount  of  hydrogen  produced  through  an 
electrolysis  process  depends  only  on  current  density,  as  stated  in 
Eq.  (5).  Since  the  implementation  of  the  compressor  power  limit 
has  no  effect  on  the  current  density  profiles,  the  variation  in 
compressor  power  is  not  expected  to  affect  the  amount  of  hydrogen 
produced. 


mH2(t) 


0. 

nF 


rt=tj 

Jt= o 


j(t)dt 


nF 


(5) 


The  total  amount  of  hydrogen  produced  over  3000  s  is  1.09  mol 
in  the  case  of  maximizing  hydrogen  production,  and  0.16  mol  in  the 
case  of  minimizing  SOEC  energy  consumption.  These  two  extreme 
cases  provide  useful  quantitative  information  on  the  range  of 
possible  behaviours  of  the  system.  The  objective  of  maximizing 
hydrogen  production  gives  a  hydrogen  production  seven  times 
larger  than  when  minimizing  energy  consumption. 

For  both  optimal  control  strategies,  the  total  energy  consump¬ 
tion  of  the  SOEC  system  (which  includes  the  SOEC  stack  and  the 
compressor)  per  unit  hydrogen  production  increases  linearly  with 
increasing  compressor  power  limit,  Pcpmax-  The  increase  in  the  total 
energy  consumption  is  mainly  due  to  the  increase  in  the 
compressor  energy  consumption  because  a  higher  Pcpmax  allows  for 
a  higher  air  ratio  and  thus  a  higher  Pep  according  to  Eqs.  (1 )  and  (2). 
The  energy  consumption  of  the  SOEC  stack  per  unit  of  hydrogen 
production  only  decreases  slightly  with  an  increase  in  the 
compressor  power  limit  Pcpmax-  This  is  because  the  bound  on  the 
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Fig.  7.  Optimal  profiles  of  (a)  current  density  and  SOEC  energy  consumption,  (b)  air 
structure,  in  the  case  of  minimizing  energy  consumption  of  the  SOEC  stack. 

compressor  power  does  not  have  any  effect  on  the  current  density 
profile,  but  only  has  a  minor  effect  on  cell  voltage,  with  a  higher  air 
ratio  (as  a  result  of  a  higher  Pcpmax)  giving  a  slightly  lower  cell 
voltage  and  in  turn  a  slightly  lower  SOEC  power  demand.  The  gap 
between  the  SOEC  energy  consumption  per  unit  hydrogen  pro¬ 
duction  Esoec  and  the  total  energy  consumption  per  unit  hydrogen 
production  £totai  becomes  bigger  at  a  larger  compressor  power 
limit.  For  example,  in  the  case  of  maximizing  hydrogen  production 
as  shown  in  Fig.  8(a),  when  a  strict  bound  is  imposed  on  the 
compressor  power,  at  Pcpmax  =  4  W,  Esoec  is  96%  of  £totai.  indicating 
that  the  SOEC  stack  is  the  dominant  energy  consumer;  when  a  loose 
bound  is  imposed  on  the  compressor  power,  at  Pcpmax  =  120  W, 
Esoec  is  only  42%  of  £totai.  meaning  that  compressor  is  the  major 
energy  consumer.  The  similar  trend  is  observed  with  the  case  of 
minimizing  the  SOEC  energy  consumption.  Therefore,  from  the 
energy  saving  point  of  view,  one  may  choose  to  operate  the  system 
with  a  small  compressor  power  limit. 

On  the  basis  of  total  energy  consumption  per  unit  hydrogen 
production,  the  control  strategy  of  minimizing  the  SOEC  energy 
consumption  does  not  bring  benefit,  compared  to  the  control 
strategy  of  maximizing  hydrogen  production.  This  is  clearly  shown 
in  Fig.  8(a)  and  (b).  For  example,  for  the  case  of  a  strict  compressor 
power  bound  when  Pcpmax  =  4  W,  the  total  energy  consumption  per 
unit  hydrogen  production  in  minimizing  SOEC  stack  energy  con¬ 
sumption  is  72.48  Wh  per  mole  H2,  compared  to  70.76  Wh  per  mole 
H2  when  maximizing  hydrogen  production  at  Pcpmax  =  4  W.  For  the 
case  of  a  loose  compressor  power  bound,  the  total  energy  con¬ 
sumption  per  unit  hydrogen  production  when  minimizing  SOEC 
stack  energy  consumption  at  Pcpmax  =  16  W  is  134.64  Wh  per  mole 


(c)  temperature  difference  and  (d)  maximum  local  temperature  gradient  in  the  solid 


H2,  compared  to  81.56  Wh  per  mole  H2  when  maximizing  hydrogen 
production  at  Pcpmax  =  16  W.  However,  when  these  numbers  are 
multiplied  by  the  amount  of  hydrogen  produced,  the  total  energy 
consumption  at  Pcpmax  =  16  W  is  20.93  Wh  for  minimum  SOEC 
stack  energy  consumption,  and  88.7  Wh  for  maximum  hydrogen 
production,  showing  a  much  larger  energy  consumption  in  the 
latter  case. 

The  choice  of  control  strategy  is  thus  a  matter  of  priority, 
depending  on  whether  the  amount  of  hydrogen  produced  is  the 
main  concern,  or  the  amount  of  energy  consumed.  Such  a  choice 
will  automatically  depend  on  factors  such  as  location,  season, 
electricity  price  and  hydrogen  price,  etc.  For  example,  one  may 
want  to  operate  the  SOEC  system  by  switching  between  the  two 
control  situations;  maximizing  hydrogen  production  when  surplus 
electricity  from  intermittent  energy  sources  is  abundant  and  the 
electricity  price  is  relatively  low,  and  minimizing  SOEC  electricity 
consumption  at  times  when  little  surplus  electricity  from  inter¬ 
mittent  energy  sources  is  available  and  the  electricity  price  is  high. 

3.4.  Minimizing  energy  consumption  by  the  air  compressor 

In  this  section  we  look  at  Scenario  B  where  the  electricity  input 
to  the  SOEC  system  (as  illustrated  in  Fig.  1 )  acts  as  a  disturbance,  to 
mimic  the  variations  of  power  over  time  found  with  intermittent 
renewable  energy  sources  such  as  wind  and  solar  energy.  The 
variations  of  power  from  these  energy  sources  can  occur  at  all 
timescales:  seconds,  minutes,  hours,  days,  months,  seasons  and 
years.  Electrical  power  systems  can  be  designed  to  cope  effectively 
with  these  variations  for  optimal  utilization  of  the  power,  and  can 
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Fig.  S.  Hydrogen  production  and  energy  consumption  as  a  function  of  maximum 
compressor  power  in  the  case  of  (a)  maximizing  hydrogen  production  and  (b)  mini¬ 
mizing  electrical  energy  consumption  of  the  SOEC  stack. 


provide  relatively  smooth  power  [27],  It  should  be  noted  that  we 
are  not  trying  to  simulate  the  ‘real’  fluctuations  in  power  supply  in 
this  paper,  but  rather  to  study  the  feasibility  of  different  control 
strategies  for  the  SOEC  system  when  there  are  fluctuations  in 
electricity  supply.  For  this  purpose,  the  electricity  disturbance  in 
the  power  supply  is  arbitrarily  simulated  as  step-wise  constant 
profile  that  reaches  a  peak  and  then  declines.  The  time  horizon  is 
chosen  to  be  3000  s,  to  be  consistent  with  the  optimal  control 
studies  for  Scenario  A,  and  also  to  represent  the  variations  in  power 
within  a  one-hour  timeframe,  as  these  have  the  most  significant 
impact  on  the  system  from  the  technological  point  of  view  [28].  In 
Fig.  9,  such  a  simulated  peak  with  five  time  intervals,  t\— 15,  is 
shown;  the  power  increases  from  20  W  to  100  W,  and  then 
decrease  to  20  W,  in  equal  increments  of  40  W. 

In  order  to  understand  the  impact  of  the  input  power  distur¬ 
bances  on  the  SOEC  system,  three  power  profiles  with  different 
shapes  (named  as  cl,  c2  and  c3),  shown  in  Fig.  9,  are  simulated  and 
used  in  minimizing  the  compressor  energy  consumption.  All  three 
power  profiles  have  step  profiles  starting  at  20  W  and  reaching  a 
maximum  at  100  W.  The  first  and  the  last  intervals  are  kept  the 
same  for  all  the  three  power  profiles,  whilst  the  three  intervals  in 
the  middle  vary  in  duration,  resulting  in  power  profiles  with  nar¬ 
rower  (c2)  to  broader  (c3)  peaks. 

For  a  given  input  power  profile,  we  set  the  optimal  control 
objective  to  minimize  the  compressor  energy  consumption,  which 
in  turn  will  lead  to  the  energy  consumption  and  therefore  hydrogen 
production  of  the  SOEC  stack  being  maximized.  Fig.  10(a)— (f)  shows 
the  optimal  control  profiles  when  cl  in  Fig.  9  is  input  as  the  total 
power  of  the  SOEC  system.  The  optimal  control  profiles  with  power 
profiles  c2  and  c3  are  not  reported  here  as  they  have  very  similar 


shapes  as  those  for  cl,  shown  in  Fig.  10.  As  can  be  seen  from 
Fig.  10(a),  by  minimizing  the  compressor  energy  consumption,  the 
majority  of  the  input  power  flows  to  the  SOEC  stack.  The  power 
profile  of  the  SOEC  stack  in  turn  determines  the  current  density 
profile,  with  similar  step  changes  in  the  range  of  1430—7600  A  itT2, 
as  shown  in  Fig.  10(b).  The  air  ratio  decreases  from  the  initial  value 
of  3.6,  reaching  its  minimum  value  of  0.4  once  the  current  density 
increases  to  4700  A  uT2  at  t  =  200  s,  and  staying  constant  at  0.4  for 
most  of  the  control  period  and  then  increasing  again  when  the 
current  density  drops  from  4700  to  1600  A  m  2  at  f  =  2300  s,  in 
order  to  maintain  the  temperature  difference  and  gradients  within 
the  constraints  specified,  as  shown  in  Fig.  10(c).  The  temperature 
difference  AT  rises  gradually  from  AT  =  -87.1  K  at  the  start  to 
AT  =  -31.5  K  at  the  peak,  and  then  gradually  falls  to  AT  =  -88.6  K  at 
the  end  of  the  time  horizon,  as  shown  in  Fig.  10(d).  The  maximum 
local  temperature  gradient  in  the  solid  structure,  however,  shows 
more  fluctuation  with  time,  as  shown  in  Fig.  10(e);  the  highest 
gradient  is  found  to  be  10  K  cm  1  when  the  current  density  drops 
from  4700  to  1600  A  nrr2;  the  lowest  gradient  is  2.5  K  cur1  during 
the  interval  when  the  current  density  is  kept  constant  at 
7600  A  itT2.  The  cumulative  hydrogen  production  profile  in 
Fig.  10(f)  shows  a  non-linear  increase  with  time;  the  total  amount  of 
hydrogen  produced  over  the  entire  time  horizon  is  0.69  mol. 

The  amount  of  hydrogen  produced  and  the  energy  consumption 
over  the  entire  time  horizon  are  summarized  in  Table  5  for  all  the 
three  input  power  profiles  cl— c3.  From  Table  5  we  can  see  that  the 
lowest  hydrogen  production  and  lowest  energy  consumption  (in 
terms  of  both  total  and  per  unit  hydrogen  production)  are  found 
with  the  input  power  profile  c2;  whilst  the  highest  hydrogen  pro¬ 
duction  and  the  highest  energy  consumption  (in  terms  of  both  total 
and  per  unit  hydrogen  production)  are  found  with  the  total  input 
power  profile  c3. 

For  all  the  cl— c3  cases  investigated  here,  the  SOEC  energy 
consumption  is  around  33  kWh  per  kg  H2,  and  the  total  energy 
consumption  per  unit  hydrogen  production  is  around  35  kWh  per 
kg  H2.  To  compare,  the  reported  total  energy  consumption  of 
Hydrogenics  and  Norsk  Hydro's  alkaline  electrolyser  system,  and 
Proton's  PEM  electrolyzer  system,  including  all  the  auxiliary  units, 
are  53.84  kWh  per  kg  H2  and  70.18  kWh  per  kg  H2  respectively  [29], 
Realistically,  the  total  energy  consumption  of  the  SOEC  system  is 
going  to  be  higher  than  35  kWh  per  kg  H2.  This  is  because:  (1)  a 
considerable  amount  of  energy  is  required  to  heat  up  the  SOEC 
stack  and  run  auxiliary  units  such  as  heat  exchangers,  which  is  not 
included  in  the  SOEC  system  studied  in  this  paper;  (2)  in  reality, 
inhomogeneity  occurs  in  the  SOEC  stack,  normally  with  cells  at 
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Fig.  10.  Optimal  profiles  of  (a)  power,  (b)  current  density,  (c)  air  ratio,  (d)  temperature  difference,  (e)  energy  consumption,  and  (f)  hydrogen  production,  in  the  case  of  minimizing 
compressor  energy  consumption  for  the  total  input  power  profile  cl  (solid  line  in  Fig.  9). 


both  ends  having  higher  voltages  (thus  higher  energy  consump¬ 
tion)  than  those  in  the  middle  of  the  stack.  Nevertheless,  an  energy 
consumption  of  35  kWh  per  kg  H2  for  an  electrolyser  system  which 
includes  the  major  energy  consumption  units  (i.e.  an  SOEC  and  a 
compressor)  is  still  competitive. 

4.  Conclusions 

A  study  of  optimal  control  strategies  for  hydrogen  production 
based  on  SOEC  technology  has  been  presented  in  this  paper,  with 
the  aim  of  offering  efficient  large-scale  system  operation  when 


coupling  an  SOEC  system  with  renewable  energy  sources.  The 
system  model  used  includes  a  ID  dynamic  SOEC  stack  model  and 
an  air  compressor  model,  to  examine  hydrogen  production  in 
relation  to  energy  consumption.  Several  scenarios  and  optimal 
control  strategies  were  considered:  in  Scenario  A,  in  which  the 
current  density  is  a  control  variable  allowed  to  vary  over  the  range 
1000—7000  A  in  2 ,  the  impact  of  maximizing  hydrogen  production 
or  minimizing  SOEC  energy  consumption  was  considered.  In  Sce¬ 
nario  B,  in  which  the  input  power  is  treated  as  a  disturbance,  with  a 
step-wise  profile  in  the  range  of  20— 100  W,  the  minimization  of  the 
air  compressor  energy  consumption  was  investigated.  Temperature 
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Table  5 

Summary  of  hydrogen  production  and  energy  consumption  for  the  control  strategy 
discussed  in  Section  3.4. 


^H2,  total 

(moles) 

Esoec 

(Wh) 

Etotal 

(Wh) 

Esoec  per  kg  H2 
(kWh  kg-1  H2) 

Etotai  per  kg  H2 
(kWh  kg-1  H2) 

cl 

0.69 

45.02 

47.78 

32.62 

34.62 

c2 

0.62 

40.71 

43.33 

32.83 

34.94 

c3 

0.81 

53.52 

56.67 

33.04 

34.98 

control  of  the  SOEC  stack  was  imposed  via  path  constraints  for  the 
overall  temperature  difference  across  the  cell  and  the  local  tem¬ 
perature  gradient,  to  link  material  properties  with  system  perfor¬ 
mance.  The  optimal  control  profiles  showed  that  these  control 
strategies  are  feasible  under  a  range  of  conditions.  It  was  shown 
that  it  is  important  to  impose  a  limit  on  the  air  compressor  power  in 
order  to  constrain  the  system  energy  consumption,  when  this  is  not 
considered  in  the  objective  function. 

The  two  control  strategies  —  maximizing  hydrogen  production 
and  minimizing  SOEC  energy  consumption  —  in  Scenario  A,  were 
compared  in  terms  of  hydrogen  production  and  energy  consump¬ 
tion.  As  expected,  the  amount  of  hydrogen  produced  within  the 
control  horizon  (3000  s)  is  much  higher  for  maximizing  hydrogen 
production  than  that  for  minimizing  SOEC  energy  consumption. 
The  energy  consumption  per  unit  hydrogen  production  (Wh  per 
mole  H2)  for  the  two  control  strategies  did  not  vary  much.  This 
provides  guidance  on  what  control  strategy  to  choose  depending 
on  location,  electricity  price  and  hydrogen  demand  etc. 

in  Scenario  B,  when  variation  in  the  input  power  is  introduced, 
the  control  strategy  of  minimizing  compressor  power  was  shown  to 
have  the  expected  effect  of  maximizing  power  sent  to  the  SOEC 
stack,  thus  maximizing  hydrogen  production.  Different  power 
profiles  were  used  as  input  to  the  SOEC  system.  The  objective 
function  chosen  was  found  to  ensure  that  the  energy  consumption 
per  unit  hydrogen  production  was  very  similar  for  the  three 
different  input  power  profiles  considered.  The  energy  consumption 
of  the  SOEC  system  studied  was  found  to  be  35  kWh  per  kg  H2  in  all 
cases.  This  is  highly  competitive  with  other  methods  of  hydrogen 
production  and  provides  a  strong  incentive  for  further  develop¬ 
ment  of  SOEC  based  technologies  for  hydrogen  production. 

Acknowledgement 

Q,  Cai  is  grateful  for  the  funding  support  from  the  EPSRC  (En¬ 
gineering  and  Physical  Sciences  Research  Council  of  UK)  platform 
grant  (EACPR  PS0693)  during  this  research. 

Appendix  A.  Supplementary  data 

Supplementary  data  related  to  this  article  can  be  found  at  http:// 
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Nomenclature 

Cp.c,  cpy 1,  cp,s,  cpj  heat  capacity  of  the  cathode  and  anode  gas 
streams,  solid  structure  and  interconnect 

(J  kg^1  K-1) 

Cp.air  heat  capacity  of  air  (J  kg  1  K  1 ) 

Ch2,  Ch2o.  Co2,CN2  concentration  of  hydrogen,  steam,  oxygen  and 
nitrogen  in  the  cathode  and  anode  gas  streams 
(mol  m  3) 

Cj  concentration  of  gas  species  i  in  the  cathode  and  anode 

gas  streams  (mol  m-3) 

Cp,  Cl  concentration  of  gas  species  i  at  the  cell  inlet  and  outlet 
(mol  m-3) 


C™,  C™.  CjpB  concentration  of  hydrogen,  steam  and  oxygen  at 
the  triple  phase  boundary  (mol  m  3) 
dh,c,  dhy \  hydraulic  diameter  of  the  cathode  and  anode  gas 
channels  (m) 

Oeff, cathode.  Deff, anode  average  effective  diffusivity  coefficient  of  the 
cathode  and  anode  (m2  s-1) 

^electrode  activation  energy  for  the  exchange  current  density  of  an 
electrode  (J  mol-1) 

Esoec.  Ecpenergy  consumption  of  the  SOEC  stack  and  the  air 
compressor  (kWh) 

/°,/c°  inlet  flow  rate  of  the  anode  and  cathode  streams  (mol  s-1) 

F  Faraday's  constant  (C  mol-1) 

hc,  h/\,  hs,  hi  height  of  the  cathode  and  anode  channels,  and 

thickness  of  the  solid  structure  and  interconnect  (m) 
j(x)  local  current  density  in  the  x  direction  (A  m-2) 
jo, cathode,  jo.anode  exchange  current  density  of  the  cathode  and 
anode  (A  m-2) 

j  average  current  density  of  a  cell  (A  m-2) 

fcc,  kA  convective  heat  transfer  coefficient  between  the  solid 

parts  of  the  cell  and  the  cathode  and  anode  gas  streams 
(J  m-2  s-1  K-1) 

^electrode  pre-exponential  factor  for  the  exchange  current  density  of 
an  electrode  (Q-1  m-2) 

L  length  of  the  SOEC  stack  in  the  z  direction  (m) 

m  air  flow  rate  into  the  compressor  (kg  s-1) 

n  the  number  of  electrons  transferred  in  an 

electrochemistry  reaction,  which  is  2  for  the  steam 
electrolysis  reaction 

Nuc,  Nua  Nusselt  number  of  the  cathode  and  anode  gas  streams 
p q2  oxygen  pressure  in  the  anode  gas  stream  (Pa) 

Pair  out.  Pair  in  the  outlet  and  inlet  air  pressure  of  the  air  compressor 

(Pa) 

Psoec.  Pc p.  Ptotai  power  input  to  the  SOEC  stack  and  the  air 

compressor,  and  system  total  power  demand  (W) 
R  reaction  rate  (mol  s-1  m-2) 

Rohm  total  resistance  of  the  cell  including  electric  and  ionic 
resistances  (Om2) 

3?  gas  constant  (J  mol-1  I<-1) 

t  time  (s) 

tf  optimal  control  time  horizon  (s) 

T  temperature  (K) 

Tc,  Ta,  Ts,  T/  temperature  of  the  cathode  and  anode  gas  streams, 
solid  structure  and  interconnect  (K) 

T J?,  Tjl  temperature  of  the  cathode  gas  stream  at  the  inlet  and  the 

outlet  (K) 

uc,  ua  velocity  of  the  cathode  and  anode  gas  streams  (m  s-1) 

U  cell  potential  (V) 

U°  standard  potential  (V) 

Urev  reversible  potential  (V) 

L/rev,TpB  reversible  potential  estimated  for  the  gas  concentrations 
at  the  triple  phase  boundary  (V) 
w  width  of  a  cell  in  the  x  direction  (m) 

x  axial  coordinate  (m) 

yo2  mole  fraction  of  oxygen  in  the  anode  gas  streams 
Yh2  hydrogen  production  rate  per  (mol  s  ') 

YH2  totai  total  hydrogen  production  over  the  entire  control  time 
horizon  (mol) 

a  transfer  coefficient 

y  ratio  of  the  specific  heat  of  air  y  =  1.4 

es,  ei  emissivity  of  the  solid  structure  and  interconnect 

Pact, cathode.  Pact, anode  activation  overpotential  losses  at  the  cathode 
and  anode  (V) 

Pconc, cathode.  Pconc.anode  concentration  overpotential  losses  at  the 
cathode  and  anode  (V) 
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Vohm  Ohmic  losses  (V) 

Vcp,  Vem  compressor  efficiency  and  compressor  drive  motor 
efficiency 

Ac,  A/i,  As,  A i  thermal  conductivity  of  the  cathode  and  anode  gas 
streams,  solid  structure  and  interconnect 
(Jnr1  s-1  K”1) 

Vi  stoichiometric  coefficient  of  gas  species  i 

vo2  stoichiometric  coefficient  of  oxygen 

Pc,  Pa,  Ps,  Pi  density  of  the  cathode  and  anode  gas  streams,  solid 
structure  and  interconnect  (kg  itT3) 

Pair  density  of  air  (kg  itT3) 

a  Stefan— Boltzmann  constant  (W  rrr2  I<~4) 

^cathode,  ^anode  electric  conductivity  of  the  cathode  and  anode 
(Q_1  rrr1) 

^electrolyte  ionic  conductivity  of  the  electrolyte  (CT1  itT1) 
t cathode.  reiectroiyte.  ranode  thickness  of  the  cathode,  electrolyte  and 
anode (m) 

\)/  air  ratio 
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